Ordered phase and scaling in Z„ models and the three-state antiferromagnetic Potts 

model in three dimensions 
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Based on a Renormalization-Group picture of Zn symmet- 
ric models in three dimensions, we derive a scaling law for the 
Zn order parameter in the ordered phase. An existing Monte 
Carlo calculation on the three-state antiferromagnetic Potts 
model, which has the effective Zq symmetry, is shown to be 
consistent with the proposed scaling law. It strongly supports 
the Renormalization-Group picture that there is a single mas- 
sive ordered phase, although an apparently rotationally sym- 
metric region in the intermediate temperature was observed 
numerically. 
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I. INTRODUCTION 

The symmetry and the dimensionality are important 
factors to determine the universality class of critical phe- 
nomena. The 0(2) symmetry is the simplest among the 
continuous symmetry, and statistical models with the 
0(2) symmetry has been studied intensively. A natu- 
ral question then would be the effect of the symmetry 
breaking from the continuous 0(2) to the discrete Zn- A 
simple spin model with Z„ symmetry is the n-state clock 
model with a Hamiltonian 



H 



o), 



(1) 



where (j, k) runs over nearest neighbors, and 9j takes 
integral multiples of 27r/n. The standard XY model with 
0(2) symmetry is defined by the Hamiltonian of the same 
form; the only difference is that 9 takes continuous values. 

The Zn symmetry is fundamentally different from 0(2) 
because of its discrete nature. On the other hand, for 
large n, it is natural to expect the Zn symmetry to have 
similar effects to that of the 0(2) symmetry. Under- 
standing these two apparently contradictory aspects is 
an interesting problem. Besides the theoretical motiva- 
tion, there are some possible experimental realizations 
of the effective Z„ symmetry. For example, the stacked 
triangular antiferromagnetic Ising (STI) model with ef- 
fective Zq symmetry may correspond to materials such 
asti CsMnls. 

In two dimensions, the phase diagram of the Z„ model 
is well understood3 in the framework of the renormal- 
ization group (RG). For n > 5, there is an intermediate 
phase between the low-temperature ordered phase with 



the spontaneously broken Z„ symmetry and the high- 
temperature disordered phase. The intermediate phase is 
0(2) symmetric and corresponds to the low-temperature 
phase of the XY model. _ 

On three dimensional (3D) case, Blankschtein et al.tl 
in 1984 proposed an RG picture of the Zg models, to 
discuss the STI model. They suggested that the transi- 
tion between the ordered and disordered phases belongs 
to the (3D) XY universality class, and that the ordered 
phase reflects the symmetry breaking to Zq in a large 
enough system. It means that there is no finite region of 
rotationally symmetric phase which is similar to the or- 
dered phase of the XY model. Unfortunately, their paper 
is apparently not widely known in the related fields. It 
might be partly because their discussion was very brief 
and not quite clear. 

In fact, there has been a long-standing controversy on 
the the three-state antiferromagnetic Potts (AFP) model 
on a simple cubic lattice, defined by the Hamiltonian 



H = 



(2) 



where aj = 0,1,2 and {j,k) runs over nearest neighbor 
pairs on a simple cubic lattice. The order parameter of 
this model is not evident. However, previous studies re- 
vealed that the low-temperature ordered phase, whi' 
is called as Broken Sublattice Symmetry (BSS) phaS' 
corresponds to a spontaneous breaking of the Zq sym- 
metry. Thr^ the effective symmetry of this model may 
be regarded^ as Zq, although it is not apparent in the 
model. It is now widely accepted that there is a phase 
transition with critical eipopents characterized by the 
3D XY universality clasaJllj, at temperature Tc ^ 1.23 
(we set the Boltzmann constant fc^ = 1.) On the other 
hand, according to numerical calculations, there appears 
to be an intermediate phase below Tc and above the 
low-tempC|i'a.ture phase. While there have been various 
proposalsETLl for the intermediate region, most reliable 
numerical results at present indicates that the interme- 
diate region appears to be rotationally symmetric phase 
which|jS|-S|imilar to the ordered phase of the 3D XY 
modeEjil3. However, the "transition" between the in- 
termediate region and the low-temperature phase is not 
well understood. According to the suggestion in Rcf. ||, 
the intermediate "phase" would be rather a crossover to 
the low-temperature massive phase. 

On the othepJiand, there has been a claim of an inter- 
mediate phaseEj also in the 6-state clock (6CL) model, 
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which has the manifest symmetry. In a recent de- 
tailed numerical study, Miyashitacd found that the inter- 
mediate region appears to have a rotationally symmetric 
character, as found in the AFP model. However, through 
a careful examination of the system size dependence, he 
concluded that it is just a crossover to the massive low- 
temperature phase, and that the rotationally symmetric 
XY phase does not exist in the thermodynamic limit. His 
conclusion is consistent with the suggestion in Ref. ^. 

In this article, based on the RG picture, we derive a 
scaling law of an order parameter which measures the ef- 
fect of symmetry breaking from 0(2) to Z„. We demon- 
strate that the Monte Carlo results on the AFP model in 
Rcf. ^ is consistent with the scaling law, supporting the 
RG picture with a single phase transition. 



II. RENORMALIZATION-GROUP PICTURE 

Since the discussion of the RG picture in Ref. ^ was 
rather brief, it would be worthwhile to present the RG 
picture here, with some clarifications and more details. 
We also make a straightforward extension to general in- 
teger n from the n = 6 case. 

A generic symmetric model may be mapped, in the 
long-distance limit, to the following ^''-type field theory 
with the Euclidean action 



?/|$P+.g|$|4-A„($"-t-$")] (3) 



with the complex field $ and its conjugate The A„- 
term is the lowest order term in $ which breaks the 
symmetry from 0(2) to Z„. The phase transition cor- 
responds to the vanishing of (the renormalized value of) 
the parameter u. The temperature T in the Z„ statistical 
system roughly corresponds to w as u ^ T — Tc where Tc 
is the critical temperature. 

In the absence of the symmetry breaking A„, the tran- 
sition belongs to the so-called 3D XY universality class. 
Its stability under the symmetry breaking to Z„ is de- 
termined by the scaling dimension of A„ at the 3D XY 
fixed point. It may be estimated with the standard e- 
expansion method. 

The lowest order result in e can be easily obtained from 
the Operator Product Expansion (OPE) cocfficicntsEj. 
As a result, we obtain the scaling dimension of A„ in 
4 — e dimensions as 



, n nin — 1) 

w„=4-n + e 1 ^ - 

yn ^ I 2 10 



(4) 



Hn is defined so that the effective strength of the pertur- 
bation A„(/) at scale I is proportional to Z^" near the XY 
fixed point. The case n = 4 is actually the special case 
N ~ P-oi the "cubic anisotropy" on the 3D 0{N) fixed 
pointE£l. Extrapolating the 0{e) result to 3D (e = 1), we 
see that the Z„ perturbation is irrelevant at the 3D XY 
fixed point for n > ric- The threshold ric is estimated 



to be 4 in 0(e) In fact, n = 2 and n = 3 corresponds 
to the 3D Ising and 3-state (ferromagnetic) Potts model, 
which do not belong to XY universality class. Thus ric 
is expected to be at least 4. This is consistent with the 
above result from 0{e). However, extrapolating the low- 
est order result in e to 3D (e = 1) is not quite reliable; 
the true value of ric might be larger than 4. On the other 
hand, we can make following observation. For n > 6, A„ 
is marginal or irrelevant at the 3D Gaussian fixed point 
(5 = 0). Thus it is natural to expect them to be irrelevant 
at the more stable 3D XY fixed point, namely ric < 6. In 
fact, the numerical observation of the 3D XY universality 
class in 6CL and AFP model strongly suggests that Ag 
is irrelevant at the XY fixed point and hence ric < 6. In 
the following, we restrict the discussion to the irrelevant 
case n > ric- 

For the 0(2) symmetric case A„ = 0, low-temperature 
phase w < is renormalized to the low-temperature fixed 
point. It describes the massless Nambu-Goldstone (NG) 
modes on the groundstate with the spontaneously broken 
0(2) symmetry. Let us call the low-temperature fixed 
point as NG fixed point. In terms of the field theory, it is 
described by the 0(2) sigma model (free massless boson 
field) 



S = 



(5) 



where (j> is the angular variable $ ~ j^je*"^. Namely, only 
the angular mode cj) remains gapless as a NG boson. In 
three dimensions, the coupling constant K renormalizes 
proportional to the scale I, and goes to infinity in the low- 
energy limit. The coupling constant may be absorbed by 
using the rescaled field 9 — \/X(0— </)o) so that the action 
is always written as / d'^x (9^0)^/2. 

Now let us consider effects of the symmetry break- 
ing A„. The symmetry breaking term can be written 
as -A„(<i>" + l>") = -A„|<i>|" cosni/). Using the rescaled 
field 9, the total effective action at scale / becomes 



5*== / d^x-{df,9f - XnK^ I d^xcoa[n{(l)o + 



e 



K 



(6) 



where the factor ^ P comes from the scale trans- 
formation of the integration measure. In the thermody- 
namic limit, we should take K —^ 00 limit. Physically, it 
means that the 0(2) symmetry is spontaneously broken 
so that the angle is fixed to some value (po in a single 
infinite system. Then the Taylor expansion of the cosine 
in 9/VK becomes valid: 

9 °° 

cos [n(</.o + -^)] = c,K'~'^'0', (7) 



3=0 



where 



2k 



C2k = {~lf^^cosn(j)o, 
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C2k+1 = ~(~^)'' (2fc+ 1)! ^^"""^O' ^^'> 

for a nonnegative integer k. The five terms j = 1, . . . 5 are 
relevant perturbations. For any value of (/>o, some of the 
coefficients Cj of these relevant terms are non- vanishing. 
We therefore conclude that, unlike the 2D case, the Zn 
perturbation is always relevant for any value of n at the 
NG fixed point. We emphasize that this conclusion is 
universal in three dimensions and independent of the mi- 
croscopic model. Shortly speaking, the Zn perturbation 
gives mass to the pseudo NG boson 9, which would be 
massless NG boson in the absence of the perturbation. In 
contrast, in two dimensions the coupling constant K of 
the free boson field theory is dimensionless, and the above 
argument does not apply. It is related to the absence of 
a spontaneous breaking of a continuous symmetry. 

We now have a global picture of the RG flow as shown 
in Fig. |l|. The phase transition between the ordered phase 
and the disordered phase is governed by the XY fixed 
point. This means that the critical exponents are iden- 
tical to those of the XY model. This is consistent with 
the numerical results. In the disordered phase above Tc, 
there will be no essential effect of the Z„ perturbation. 
However, the nature of the ordered phase is more inter- 
esting. The Zn perturbation A„ is eventually enhanced 
in the ordered phase below Tc. It means that all region 
below Tc belong to the massive phase with the sponta- 
neously broken Z„ symmetry. There is no rotationally 
symmetric intermediate phase, unlike the 2D case. Only 
a precisely 0(2) symmetric model with A„ = is renor- 
malized to the NG fixed point below Tc, corresponding 
to the rotationally symmetric low-temperature phase. 

An interesting aspect of the RG flow diagram is that 
the Zn perturbation is irrelevant at the 3D XY fixed point 
but is relevant at the low-temperature NG fixed point. 
This could be related to a nontrivial system size depen- 
dence found, in a Monte Carlo Renormalization Group 
calculationll3. For T slightly less than Tc, the symme- 
try breaking perturbation A„ is renormalized to a small 
value by the RG flow, and remains small until the RG 
flow reaches near the NG fixed point. It means that the 
mass of the pseudo-NG bosons is suppressed by the fluc- 
tuation effect. At a flnite scale (for example in a flnite 
size system), the ordered phase near Tc is very similar 
to the low-temperature phase of the XY model. This 
naturally explains the numerical observation of the ap- 
parently rotationally symmetric "phase" in 6CL or the 
AFP model. For larger n, the mass is more suppressed, 
and the low-temperature side of the transition appears 
to be O (2) symmetric until the system size becomes very 
large. However, for any flnite n, the low-temperature 
side of the transition T < Tc is not truly massless nor 
0(2) symmetric in the thermodynamic limit, as already 
pointed out. 



III. SCALING LAW IN THE ORDERED PHASE 

Based on the RG picture, we derive a scaling law on 
an order parameter On which characterizes the symme- 
try breaking from the 0(2) to Z„ symmetry. There are 
various passible deflnitions of On- On the 6CL model, 
Miyashitao numerically measured an order parameter A 
which corresponds to the effective barrier height. On the 
AFP model, Heilmann, Wang and SwendsenU studied 
{4>6), which is the Fourier transform of the angle distribu- 
tion density of average spins. The following consideration 
apply to the both cases. 

For large enough L and T slightly lower than Tc we 
divide the RG flow to three stages, as shown in Fig. || 

(i) The RG flow near the 3D XY flxed point. The sym- 

metry breaking A„ is irrelevant, and is renormalized 
proportional to l~\y^\ at length scale /. 

(ii) The RG flow from the neighborhood of the 3D XY 

fixed point to the NG fixed point. For simplicity, 
we assume that the symmetry breaking A„ is un- 
changed in this stage. 

(iii) The RG flow near the NG flxed point. A„ is rele- 
vant, giving a mass to the NG boson. 

The length scale Ic, at which the crossover from Stage 
(i) to (ii) occurs, is given by Ic ^ const. (Tc — T)^" , 
where v is the correlation length exponent of the 3D XY 
universality class. Thus, at the crossover, 

A„ -const.(Tc-r)''ly"l. (9) 

This also gives the effective value of the perturbation A„ 
at the crossover from Stage (ii) to Stage (iii). 

In the presence of the Z„ perturbation, the spin config- 
uration would be dominated by the ordered regions which 
are separated by domain walls in a large system. The 
free energy costed by the domain walls is proportional 
to their area, which scales as for the system size jL- 
Therefore the effective "barrier height" is proportionalta 
to L^. Combining this with eq. (^, we conclude that the 
order parameter is a function of a single scaling variable: 

On^ I{cL\Tc-TY\y-\), (10) 

where c is a constant. The function / is universal, but 
of course depends on the definition of the On- While 
the scaling by was used in Ref. we find that the 
temperature dependence of the order parameter is also 
governed by a scaling. Interestingly, the exponent v\yn\ 
is completely determined by the 3D XY fixed point. 



IV. COMPARISON WITH THE NUMERICAL 
RESULTS 

In the numerical study0 of the AFP model, they 
claimed the existence of the intermediate phase, in which 
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the order parameter (c^e) is very small even for relatively 
large lattice (upto L = 64). However, we re-analyze 
their data to demonstrate the scaling relation (|0|), and 
hence the validity oL the RG picture. In FigT^, we 
show the data takcnllll from Fig. 3 of Rcf. 111]. We chose 
v\yQ\ = 4.8 to give the best scaling. The data for vari- 
ous temperature and various system sizes fall remarkably 
into a single curve as a function of the scaling variable 
X = cL^(rc — T)"^!^"!. This supports the proposed scahng 
relation (p^. Furthermore, if we approximate the effec- 
tive potential by — a;cos60, the scaling function is given 
by 

_ cos (60)6^^°'- (6^) ^ h(x)_ 

•^^""^ /d0e--''(6'^) /o(x)' ^ ' 

where /„ is the modified Bessel function. Choosing c ~ 
0.025, the scaled data agree with this simple function 
rather well. We note that the data appear to deviate 
from the scaling law for small x. This may be due to the 
insufficient system size L or the relatively large statistical 
error. 

We emphasize that the present scaling relation is a 
strong evidence of the single phase transition at the tem- 
perature Tc- In contrast, the scaling of the "spontaneous 
magnetization" p = |$| oc {Tc — T)^ does not distin- 
guish our picture and the "intermediate phase" scenario 
of Ref. 

On the other hand, the scaling function f{x) for A in 
Ref. ^ is linear in x by definition. He indeed found that 
A is scaled by L^. However he did not discuss the tem- 
perature depedence. We have attempted to analyze the 
data in Figs. 6 and 7 in Ref. to find that they are 
roughly consistent with our scaling law ( [l0| ) with the ex- 
ponent vluel ~ 4. The estimate is difficult because there 
are only small number of temperature points available in 
Ref. nj. According to our picture, the exponent J^lj/el is a 
universal quantity determined by the 3D XY universality 
class. Considering the available data, the above results 
on AFP and 6CL models are consistent with the uni- 
versality hypothesis, although not conclusive. It would 
be interesting to obtain more numerical data on these 
models to check our scaling law (|lO|). 

The exponent v has been determined as ~ 0.67 for 
the 3D XY universality classta. Combining with the 
above estimates of vluel, Ivel is estimated as about 6. 
Unfortunately, the irrelevant eigenvalue ye has not been 
much discussed in the literature. The lowest order re- 
sult (^) in the e-expansion gives \ye\ = 3, which is not 
quite consistent with the numerical estimate. However, 
it is perhaps not surprising to obtain an inaccurate re- 
sult in the lowest order of the e-expansion. It would be 
interesting to carry out the calculation to higher orders 
in e, or to estimate yg by other means. 



V. CONCLUSION AND DISCUSSIONS 

In this article, we clarified a RG picture of phase struc- 
ture of 3D Zn symmetric models, which was introduced 
earlier by Blankschtein et al.0. There is no finite region of 
intermediate phase with a (spontaneously broken) 0(2) 
symmetry, but only a crossover to a massive phase where 
the discreteness of Z„ is relevant. Based on the RG pic- 
ture, we have derived a scaling law of the order paraipetcr 
in the 3D Zn models The existing Monte Carlo datallJ on 
the AFP model, which was used to claim the intermedi- 
ate phase, was shown to be consistent with the scaling 
law. Thus we conclude that the RG picture is valid on the 
AFP model, and there is only one transition at Tc ~ 1.23 
with the 3D XY universality class. 

We would like to make a few final remarks. Firstly, we 
note that the RG argument used in the present article 
does not contradict to the transition of other than XY 
universality class, because only the local stability of the 
XY fixed point was discussed. It is possible that a lattice 
model with Z„ symmetry is renormalized to another (un- 
known) RG fixed point. Actually, it appears somewhat 
controversial whether the transitioip_n£ the STI model 
belong to the XY universality classEfro. On the other 
hand, the available numerical results strongly supports 
that the 6CL and AFP models at the critical temper- 
ature are renormalized into the XY fixed point. Once 
the transition is known to be XY universality class, the 
RG picture and the scaling law discussed in this article 
should apply to the ordered phase, for the temperature 
slightly below the critical point. 

Secondly, as discussed in Ref. ^, the "bare" value of A„ 
(at a small length scale) may have opposite sign in some 
circumstances. Namely, the minima and maxima of the 
potential of are swapped. In such a case, the ordered 
phase may correspond to the Permutationally Symmet- 
ric Sublattice (PSS) phase proposed in Ref. || for the 
AFP model, or the Incompletely Ordered Phase (lOP) 
proposed in Ref. ^ for the 6CL model. In the vicinity 
of the critical point, the temperature dependence of the 
bare A„ is not essential because the leading dependence 
on the temperature is determined by the critical effect, 
as shown in eq. (|lO|). However, it may be important in 
a wider temperature range. In particular, if the bare A„ 
changes sign at some temperature lower than Tc, we 
would have a transitionH at T^. Such a transition would 
be controlled by the NG fixed point. The existing numer- 
ical data indicates that there is no such phase transition 
in the standard AFP model on the simple cubic lattice or 
in the standard 6CL model. However, such .a .transition 
may be possible iijt,some modified modelsEiE3. In fact, 
Blankschtein et ala argued it to exist in the STI model. 
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NG fixed point 3D XY fixed point 

FIG. 1. The RG fiow diagram of the Zn models, projected 
onto the two-dimensional parameter space spanned by u and 
An. The Zn perturbation A„ is irrelevant at the 3D XY fixed 
point, but is relevant at the NG fixed point. For T slightly less 
than Tc, the RG flow is divided into the three stages (i),(ii) 
and (iii). 
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FIG. 2. The order parameter ((pe) taken from Ref.t-!l. They 
are scaled by a; = cL^{Tc — T)"'^"', for various system sizes 
and temperatures. The data are consistent with the scaling 
law ( p^ with the exponent i^lj/el ~ 4.8. They also agree 
with the approximate scaling function f{x) = Ii(x)/Io{x), 
for c = 0.025. 
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